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The simulation model of a change of O, 0 2 and 0 3 oxygen component content at the 
heights 50 - 140 km under electron precipitation influence is presented here. Parameters of the air 
mass vertical transfer are introduced into the model. Calculation results showed that after the 
intense precipitation in ~ 3 days the ozone content increased at heights ~ 80 km. These results agree 
with the experimentally found effects of the content change under precipitation conditions. 

Introduction 

The experiments during solar proton events (HEATH, 1977; LIPPER, 1985; THOMAS, 
1983) show a decrease of ozone content at the flare initial phase. A further behaviour of ozone 
concentration in the stratosphere and mesosphere differ from each other. For example, a total 
ozone number at height ~ 37 km after the August 1972 proton flare (HEATH et al„ 1977) did not 
recover for several weeks. In the mesosphere the ozone concentration recovery up to its initial 
values occurred within a short time and then the further growth of [0 3 ] in some events was 
observed. During the solar proton events in July and September 1974 (LIPPER et al., 1985) at 80 
km the difference between the proton flux maximum and the ozone density minimum is ~ 36 hours 
and the maximum ozone content is observed in 4 - 6 days after the flare maximum. For the first 
time a possibility of ozone concentration increase at heights between 70 and 85 km was 
theoretically predicted by CRUTZEN (1980) by an example of the solar flare in August 1972. The 
ozone content increase after the intensive precipitation of high energy electrons on March 23, 1974 
(SKRYABIN et al., 1977) is shown. In this case the electron precipitation with average energy 
250 keV occurred as three splashes each of them lasted 15-30 min. At this time at 120 - 200 km 
height there was a transfer intensification downwards of atmospheric masses. A maximum of 
ozone content increase in the atmosphere column occurred 1 - 3 days after the precipitations and it 
was ~ 7 - 9% of the initial value of [0 3 ] (SKRYABIN et al., 1977). 

Calculation Method 

To find a "pure" contribution of precipitating electrons into the mesospheric ozone content, 
the polar night conditions are considered, i.e., the time when the short wave radiation of the Sun is 
absent. It is interesting to determine a maximum increase of ozone content in the mesosphere 
under the influence of precipitating electrons. 

It is known that under precipitation conditions in the mesosphere O atoms are generated. 
Thus, it is shown (SOCHNEV, 1977) that their formation velocity is determined by the total 
velocity of ionoformation q(t) and by concentration of the atmosphere components at given height 
(H): 

, x 3.22 [N 2 ] H[0 2 ] + 1.12 [0] 

Q ( t) 

1.15 [ N*] + 1.5 [0*3+ 0.56 [0l 

The electron flux with the exponential energy spectrum is considered here: 

J(e) = A( t) • exp (- e/?o), 
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where A(t) is a coefficient which takes into account a temporal change of precipitating electron flux 
magnitude and reaches 10^ in maximum, t is time in seconds. It is natural that such huge fluxes 
cannot exist for a long time but they can occur repeatedly within a day. As a result, such 
precipitation splashes in the auroral zone can last in total several tens of minutes. Let us assume 
that 

A(t)=B-t'exp(-t a / c ), ^ 

where B and C are constants being 5.5* 10 6 and 180000, respectively. At these values A(t) reaches 
its maximum value (10 9 cm” 2 s _1 keV -1 ) in 5 min. The average precipitation time in this case will 
be - 6 min. 

Atomic oxygen disappears in general in the following reactions (McEWAN and 
PHILLIPS, 1975; CHAMBERLAIN, 1978): 

0 + 0 M - 0 2 + M, ^*3.6 • 10 * • L. exp(- (4) 


5 1 U \ 6/r 

OfO^M - OjtM, k 2 - 1 , 1 • 10 -exp^— J,cm/S 


o+o 


, . ,.-11 ( 2150 

K 3 - 1.1-10 • exp [- ~r~ 


V s . (6) 


here M is any third particle, Kj, K2, K3 are velocity constants of the appropriate reactions. 

The reactions (4-6) take place in the atmosphere continuously. Hence, there are the 
recovery processes of [O], [OJ and [O3] which are taken to be constant at the given height and 
they are equal to the change velocity of these air components obtained from (4-6) under quiet 
conditions. These velocities are denoted by D(O), D(C>2) and D(C>3). 

The equations describing the distribution of oxygen components with the account of the 
transfer processes are as follows (SOSIN, 1985): 

3 [a0 ] 2 

=Qlt ^D(0) “2 k,, [01 [M] - [0][0 2 ] [M]~ K 5 [0][Q 3 1 +* 


a [. % 3U0] / D + K dr D K \ r n1 ] 

+ ™(d+k) — — + +■ T- + — +‘vr)U0H (7) 

3z l V T h 0 H av ) j , 


3[a0o] 


gj. - 1)(0 2 ) — i-2k 5 lo![o 5 1 +-k 1 Lo 1 [Ml-K 2 [0]t0 2 lLM]-b 


9 J . 3[*0 2 1 (T> + K dT D K „ r \ r n n 1 



344 


- D(0 3 ) t K a [0][0 2 ][M] ' k 3 [0H0 5 1 + 
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here [AO], [A0 2 ] and [A0 3 ] are increments of O, 0 2 , and 0 3 concentrations with respect to their 
initial contents before the precipitation, T is the temperature by an absolute scale in winter 
(KASTING, 1981), H 0 , H 02 , Hq 3 , and H av are the homogeneous atmosphere heights for O, 0 2 , 
0 3 and for the component with the average molecular mass, respectively. To find the average 
molecular mass the atmosphere was considered to consist of N 2 , 0 2 , 0 3 and O molecules. The 

concentration values at heights > 80 km are taken according to (KASTING, 1981). At heights < 
80 km the night concentrations of N 2 , 0 2> and 0 3 were used on(McEWAN and PHILLIPS, 1975) 
and of O on (HUNT, 1973). The transfer coefficients D, K, and V were taken from 
CHAMBERLAIN (1978); LETTAN (1951); McEWAN and PHILLIPS (1975). The calculations 
were carried out by the method of probabilistic transitions (SKRYABIN, 1985) for the transfer 
parameters 0.5 K; 0.5 D; 0.5 V and 2 K, 2 D, 2 V. 

Discussion 

In Figure 1 the simulation results of [AO], [A0 2 ] and [A0 3 ] at various times after the 
electron precipitation in a form of 6-min splash are presented. From the figure it is seen that when 
the transfer parameters increase causing the intensification of the air mass displacement 

downwards, the [AO] decreases and the [A0 3 ] increases. In Figure 2 the changes of [AO] and 
[A0 3 ] in the atmosphere column are presented. With the change of transfer parameters by 0.5, 
once and twice the maximum values of [A0 3 ] are 2*10 15 , 5*10 15 , and 9*10i 5 cm- 2 , respectively. 
It equals the change of ozone content from a quiet level by 0.3, 0.76, and 1.4%. For 30-min 
precipitation they will be 1.5, 3.8, and 7%. And the higher is the arrival velocity of the air masses 

into the lower layer, the earlier is reached the content maximum of [A0 3 ]. As is seen, such 
changes of ozone content by the order of magnitude correspond to the experimentally obtained 
effects (SKRYABIN et al., 1977). 
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Figure 1 Distributions of [AO], [AO 2 ] and [AO 3 ] 
on heights at various times: 

30 min 

6 hours 

. . 1 day 

3 days 

x x 6 days 

total content of ozone in the 

upper atmosphere on CIRA 
72 model. 

a) Coefficients of the turbulent (K), molecular 
(D) diffusions and the velocity of the averaged 
vertical wind in winter (V) according to 
CHAMBERLAIN, 1978; LETTAN, 1951; 
McEWAN and PHILLIPS, 1975);. 

(b) K, D and V are twice increased; 

(c) K, D and V are twice decreased. 
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Figure 2. The content of [AO] and [AO 3 ] in 
the atmosphere column a single cross section 
versus time. Figures on curves show K, D 
and V increase. 
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